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ABSTRACT 

To study the influence on overall deformations of the time-dependent constitutive properties 
of fiber-reinforced polymeric matrix composite materials being considered for use in orbiting 
precision segmented reflectors, simple sandwich beam models are developed. The beam 
models include layers representing the face sheets, the core, and the adhesive bonding of 
the face sheets to the core. A three-layer model lumps the adhesive layers with the face 
sheets or core, while a five-layer model considers the adhesive layers explicitly. The defor- 
mation response of the three-layer and five-layer sandwich beam models to a midspan point 
load is studied. This elementary loading leads to a simple analysis, and it is easy to create 
this loading in the laboratory. Using the correspondence principle of viscoelasticity, the 
models representing the elastic behavior of the two beams are transformed into time- 
dependent models. Representative cases of time-dependent material behavior for the face- 
sheet material, the core material, and the adhesive are used to evaluate the influence of 
these constituents being time-dependent on the deformations of the beam. As an example 
of the results presented, if it assumed that, as a worst case, the polymer-dominated shear 
properties of the core behave as a Maxwell fluid such that under constant shear stress the 
shear strain increases by a factor of 10 in 20 years, then it is shown that the beam deflection 
increases by a factor of 1.4 during that time. In addition to quantitative conclusions, several 
assumptions are discussed which simplify the analyses for use with more complicated ma- 
terial models. Finally, it is shown that the simpler three-layer model suffices in many situ- 
ations. 


INTRODUCTION 

As part of the development phase of the use of composite materials for long-duration space 
applications, it would be useful to have a simple analytical tool which models the important 
features of sandwich construction and allows for the evaluation of the influence of the time- 
dependent behavior of the various constituents in the construction; namely the face sheets, 
the core, and the adhesive. In addition, it would be beneficial to have a simple laboratory 
specimen which could be used to gather empirical data regarding time-dependent material 
behavior and screen candidate materials. Both of these requirements can be satisfied to a 
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large degree by considering sandwich beams, and by studying their time-dependent behav 
ior as a function of the time-dependent behavior of the constituents. Beams are one 
dimensional in nature, leading to somewhat simpler analyses than plate or shell specimens 
and laboratory loading of beams is generally direct and free of unwanted secondary effects 
Three- or four-point loadings are examples of simple yet effective loadings. This paper dis- 
cusses the time-dependent response of both three- and five-layer symmetric sandwich 
beams. The five-layer model includes the face sheets and the core, and the adhesive bond- 
ing these constituents together. The three-layer model lumps the adhesive layers into either 
the face sheets or the core. The response of these two models to a three-point loading is 
considered. Numerical predictions regarding the deflections over a 20 year time span are 
made in the context of the various constituents of the sandwich construction exhibiting 
time-dependent behavior. Since 20 year data are not available, the behavior of the constit 
uents is hypothesized. For more complete details of the work discussed, the reader should 
consult ref. 1 . 


DEVELOPMENT OF THE THREE-LAYER MODEL 
Nomenclature and Problem Definition 

In fig. 1 the three-layer model is described, as is the loading. The beams considered are oi 
length L. The three-point bending load consists of a simply supported beam with a point load 
P at midspan. Because of symmetry about the midspan, this loading is studied here as e 
cantilever beam with a tip load P/2, the center of the simply-supported beam being clamped 
in the analogous cantilever problem. The extensional moduli in the x direction of the face 
sheet and core are denoted as E, and E 2 , respectively. The shear moduli in the x — z plane 
are denoted as G, and G 2 . The thickness of the face sheets is t, and that of the core 2h. The 
overall thickness is 2H. 

Equations Governing Elastic Response 

For this problem it is assumed that the stress components cr y , a z , t^, and r yz are zero 
Hence the elastic stress-strain relation is given by 

(7 X =E ^ T XZ = G Vx7 • 0 .' 

In the above, for a particular layer, E is the extensional modulus in the x direction and G the 
shear modulus in the x — z plane. 

The assumed displacement field for the three-layer beam is also illustrated in fig. 1. In par- 
ticular, the sandwich cross section is assumed to displace uniformly as-a-whole in the x di- 
rection an amount u°(x) and downward as-a-whole an amount w°(x). In addition, the cross 
sections of the face sheet layers are assumed to rotate independently of the cross sectior 
of the core; the angles of rotation, <// and </>, respectively, being defined in the figure. With 
these kinematics, the displacement field is written as 

u°(x) -F h </>(x) — (z + h)i//(x) ( H < z < -h) 
u(x,z) = u°(x) - z</>(x) ( -h < z < +h) 

u°(x) - htf>(x) - (z - h)t//(x) ( Fh < z < +H) * 2 ' 

w(x,z) = w°(x) , all z . 

Hereafter, for convenience, the superscript o will be dropped and hence the strains required 
for use in the stress-strain relation are given by 
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To be determined are u(x), </>(x), i//(x), and w(x). For the viscoelastic problem the functions 
are time-dependent and should be written as u(x,t) 


The total potential energy used for determining the elastic response simplifies here to 


n = t r 1 2 « ^ + ° y ^ dxdx ■ t w ( t ) • <4) 

~~2 

where the limits on the integral reflect the fact that the analog cantilever problem is being 
considered, the tip being loaded with known load P/2. Substituting the moduli and strains 
into eq. 4 and integrating with respect to z leads to 



The constants Ci 



are as follows: 


c o — + ^2^ 

c, = h 2 (E,t, -I yE 2 h) 

c 3 = E,hl 2 

1 c *3 

c 6 3 ^1^1 
c 7 = 2G 1 t 1 
Cg 2G 2 h . 


(5) 


( 6 ) 


From eq. 5, by taking the first variation and integrating by parts, the governing equations can 
be shown to be 
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The boundary conditions associated with the variational statement are 
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at x = 0 
u = 0 

0 = 0 
i// = 0 

w = 0 


at x : 


2Cn ~7~ = 0 


2c 


0 dx 
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(c 7 + c 9 ) 


dw 
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c 7 > /' - C 90= -5- . 


On the right side, the 1st and 4th terms can be identified with inplane and shear force re- 
sultants, respectively, while the 2nd and 3rd terms are moment resultants. 

The equation for u(x) is decoupled from the other three equations. Here attention will bo 
focused on the displacement w(x) and thus the first equation will not be considered further. 
The solutions for the other three displacement variables are taken to be of the form 


Ay 


0(x) = fe 
0(x) = se 
w(x) = we 


A* 


Ax 


(91 


Substituting these forms into the last three of eq. 7 results in the following polynomial that 
must be satisfied by A: 

>i 6 (4Ci c g — c 3 )(c 7 -f- c 9 ) -- A (2c 7 c 9 )(c, t c 3 + c 6 ) = 0 . (10) 

The roots to this equation, and application of the boundary conditions, lead to 

O ry 

w(x) = w 3 x + w 2 x 4- w,x + w 0 4- w 5 sinhUx) I w 6 cosh(/tx) 

2 6(2c, 4- c 3 ) 

0(x) = 3w 3 x 4- 2w 2 x 4- [w, 4- - w 3 ] ( w 5 A 6 cosh(^x) 4- w 6 A 6 sinh(Ax) 


2 6(c 3 4- 2c 6 ) 

0(x) = 3w 3 x 4- 2w 2 x 4- [wj 4- w 3 ) 


w 5 B 6 cosh(yix) 4 w 6 B 6 sinh(^x) 


(11 


where 
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w 0 = w 5 

p (* 4 

W ’ = 2(R, 4 R 2 f V C9 

PL 

W2 " 8(R, + R 2 ) 

-P 

W 3 " 1 2(R ! + R 2 ) 

-P(c 7 R! - c 9 R 2 ) 2 

Wc = — 2 

2c 7 c 9 A(c 7 4- c 9 )(R, P R 2 ) 
w 6 = - w 5 , 

with 

20709 ( 0 , -+• C 3 P c g) 

(40, c 6 ~ c 3)( c 7 F ^ 9 ) (13) 

R, = 2c 1 -P c 3 
R 2 = c 3 + ^ C 6 • 




Interest here will focus on the response at the tip of the cantilever, this being representative 
of what the beam is doing. In that regard, using x = L/2, eqs. 11 and 12, 

/ L \ _ PL 3 , PL ( 

w ti P W V 2 ) 48(R, + R 2 ) 4(R, + R 2 ) 2 V Cs 

P(c 7 R 1 - c 9 R 2 ) 2 / 

2c 7 c g A(c 7 + c 9 )(R 1 + R 2 ) 2 

It has been found that for a very wide range of values of elastic properties the last term 
contributes very little. It is thus dropped as it considerably simplifies the algebra. 



Time-Dependent Behavior of the Three-Layer Beam 


In the present study, for obtaining an understanding of the overall effects of \he tirn 
dependence of the various constituents, and at the same time considering a worst-case sce- 
nario, the constituents are modelled as Maxwell fluids. With a Maxwell fluid for a constant 
level of applied stress the material strains indefinitely. For a Maxwell fluid the constitutive 

equation takes the form 


(7 


+ p,(7 = q,j: 


(15) 


where it is understood that a can represent a normal or a shear stress and e can represent 
an extensional or a shear strain. Taking the Laplace transform of both sides results in 


Pts')ff(s') = Q (s)r(s) 


MfVl 


where <r(s) and S(s) are the Laplace transforms of the stress and strain functions of time, re 
spectively. In the above 
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( 17 ) 


P(s) = 1 + 

Q(s) = q t s . 

It should be noted that the subscript 1 on p, and q ( in eqs. 15, 17, etc. have nothing to do with 
the face sheets. The nomenclature p, and q, are standard for viscoelastic materials. See ret. 
2, for example. 

If a particular constituent is considered to be time-dependent, then the application of the 
correspondence principle for the sandwich beam problem at hand calls for replacing, in the 
formulation for the static elastic response, the elastic modulus of the constituent with the ra- 
tio <Q>(s)/P(s), i.e., 

E or G-*Q(s)/P(s) . (18) 

The static load P is then assumed to be applied at time zero in a stepwise fashion so the load 
must be replaced with its Laplace transform, i.e., 

P-^P/s (19) 

The resulting expression is then the Laplace transform of the time-dependent response of the 
beam. Taking the inverse transform yields the response of the beam as a function of time. 

To estimate the long-term effect of time-dependent behavior in the sandwich beams here, it 
will be assumed that if a constituent exhibits time-dependent behavior, the strain of the con- 
stituent will increase by a certain factor in 20 years if subjected to a constant stress. Two 
levels of time dependency will be studied, one considered worst-case time dependency and 
the other considered minimal time dependency. The displacement response of the tip for the 
20 year period for these two levels of time dependency will then be computed. 

As an example, consider the following: With the face sheet material properties being con 
trolled to a large extent by fiber properties, a large degree of time-dependent behavior is 
unexpected. Hence, it will be assumed that the face sheet material, in the extreme case, 
exhibits a 10% increase in extensional strain when subjected to a constant stress for 20 
years. As a result, for the face sheet material 

q,=200E 1 and p 1 =200 years . (20) 

Using eq. 17, the substitutions indicated by eqs. 18 and 19 can now take place in eq. 14 (with 
the last term dropped). Performing the inverse transformation results in an expression for 
the time-dependent tip deflection, namely, 



where 
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( 22 ) 


F = t^ht, + J-t?) 

G = Aq 1 + Bp 1 
H = Cq, + Dp! 

and W(j P represents the elastic response as given in eq. 14. Consider a sandwich beam with 
quarts epoxy face sheets in an 8 layer quasi-isotropic lay-up and a honeycomb core. Table 
1 illustrates nominal constituent elastic properties. Note that table 1 inc udes mformation i for 
the five-layer model to be discussed shortly. For the three-layer model, since the adhesive 
layer is so thin, the elastic properties of the core are not adjusted to account for lumping the 
adhesive into the core. Only the thickness of the core is adjusted. 


Table 1. Nominal Material Properties and Geometries for Beam Model 


Face Sheets 

Honeycomb Core 

Adhesive 

E, = 2.5E6 psi 

E 3 = 1 .0E3 psi 

E 2 = 0.5E6 psi 

G, = 0.96E6 psi 

G 3 = 29E3 psi 

G 2 = 0.179E6 psi 

t, = 0.040 in. 

h = 0.255 in., 
3-layer model 

t 2 = 0, 

3-layer model 


h = 0.250 in., 
5-layer model 

t 2 = 0.005 in., 
5-layer model 


From eq. 21 it can be seen that the time dependency is exponential in form but, as shown in 
fiq 2 over the 20 year period it appears linear. For the quartz-epoxy/honeycomb sandwic , 
the percent increase in tip deflection, relative to the static elastic deflection at t-0, is illus- 
trated in fig. 2 for both the worst-case, or maximum, time dependency and the minimal 
time-dependency of the face-sheet material. Minimal time dependency is defined to be the 
case when the face-sheet material exhibits only a 1% increase in extensional strain when 
subjected to a constant tensile stress for 20 years. For the worst case it is seen that the 
cantilever tip deflection increases by about 8% in 20 years. For the case of minimal time 
dependency, the tip deflection increases just under 1% in 20 years. These numerical values 
reflect an almost one-to-one relationship between face-sheet material extensional properties 

and tip deflection. 
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As another example, if the core exhibits the behavior of a Maxwell fluid, the appropriate 
substitutions of eqs. 17, 18, and 19 into eq. 14 leads to an expression for the time-dependent 
tip deflection as 


w ti P 0) = w ti P + 


PL 

4(R, + R 2 ) 




J 

qi 


As can be seen, the tip deflection is linear with time. 


(23) 


Assuming for the case of maximum time dependency that the core strain increases by a 
factor of 10 in 20 years, and for the case of minimum time dependency that the core strain 
increases only by a factor of 2 in 20 years, eq. 23 leads to the results shown in fig. 3. It is 
seen that these cases lead to considerable tip deflection over a 20 year period. 


DEVELOPMENT OF THE FIVE-LAYER MODEL 


Nomenclature and Problem Definition 

To explicitly include the adhesive layer, a five-layer model is necessary. The nomenclature! 
and Kinematics of the five-layer model are shown in fig. 4. The properties of the face sheets 
are subscripted with a 1, the properties of the adhesive subscripted with a 2, and the prop 
erties of the core subscripted with a 3. Total sandwich thickness is again 2H. 

Equations Governing Elastic Response 

The sandwich cross section is again assumed to displace uniformly as-a-whole in the x di 
rection an amount u°(x) and downward an amount w°(x). The cross sections of the face shee 
layers are assumed to rotate independently of the cross sections of the adhesive layers, and 
the core has its own cross-section rotation. The displacement field is thus given by 


u(x,z) = 


u°(x) + ha(x) -l- t 2 /?(x) - (z -l- h l t 2 )y(x) 
u°(x) -I- ha(x) - (z + h)/?(x) 
u°(x) - za(x) 

u°(x) — ha(x) - (z — h)/?(x) 

u°(x) - ha(x) - t 2 /?(x) - (z - h t 2 )y(x) 


H<z< — H + t, 

- H + t,<z< — h 

— h<z<h 
h<z<H — t t 
H -t,<z<H 


w(x,z) = w°(x) . 


(24 1 


There are five kinematic quantities to be determined with the five-layer model. Dropping the 
superscript o, the strains are given by 
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(25) 


Substituting these strains into the total potential energy, eq. 4, leads to the following differ- 
ential equations for the kinematic variables: 
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c 0 = E 1*1 + E 2 { 2 + E 3 h 

Cl = h 2 (E 1 t 1 + E 2 t 2 4 yE 3 h) 

c 2 = ht 2 (2E 1 t 1 + E 2 t 2 ) 

c 3 = ^l( E 1^l) 

C 4 = + -J- Ejtjl 

c 5 = ^1^2( E 1^1 ) 

C 6 = ' 5 (yE,t,) 
c 7 = 2 G 1 t 1 
c 8 = 2G 2 t 2 
c 9 = 2G 3 h . 


The relevant boundary conditions are 


at x 

= 0 : 

at x = 

L . 
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2 c 0 
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(27 1 


(28) 


As before, the equation governing u°(x) decouples from the other four and it can be disre 
garded at this time. Solution of the equations for ot(x) , /?(x) , y(x) , and w(x) follows the 
procedure for the three-layer beam. The algebra, however, is considerably more involved. 
As with the three-layer beam, approximations are used to simplify the algebra, particularly 
for application of the correspondence principle. Through these approximations, the ex 
pression for the tip deflection of the five-layer beam takes the form 


W 


tip 



PL 3 

48(R 1 4 R 2 4 R 3 ) 


4 


PL 

4(R, 4 R 2 4 R 3 ) 2 



(29) 


where 


R 1 = 2c 1 4 c 2 4 c 3 

R 2 = c 2 4 2c 4 4 c 5 (30) 

^3 = C 3 + C 5 + 2C 6 . 


This equation is the analog of eq. 14 for the three-layer beam. 
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Time-Dependent Response of the Five-Layer Beam 


Numerical studies with the five-layer beam indicate that little new information, relative to 
what can be learned with the three-layer model, is obtained by using the five-layer model to 
study the influence of time-dependent face sheet and core properties. However, since the 
five-layer model includes the additional feature of the adhesive layer, it is of value to discuss 
the influence of time-dependent adhesive properties. Assuming the adhesive shear proper- 
ties behave as a Maxwell fluid, using the expression for the tip deflection, eq. 29, performing 
the steps given by eqs. 18 and 19, and taking the inverse transformation, the tip deflection 
as a function of time for a viscoelastic adhesive layer is 


w tip 0) = w tip + 


PL 

4(R, + R 2 + R 3 ) 2 



(31) 


This expression also depends on time in a linear manner. Using the material properties as 
given in table 1, the time-histories of the percent increase in tip deflection for the cases of the 
adhesive shear strain doubling and increasing ten-fold in 20 years are shown in fig. 5. For 
the worst-case condition, the tip deflection changes by less than 0.5% in 20 years. 


CONCLUSIONS 

Presented has been the development of models to be used in evaluating the influence of 
time-dependent face sheet, core, and adhesive constitutive properties on the overall defor- 
mations of sandwich beams. The study has its origins in the need to understand the time- 
dependent deformations of orbiting precision segmented reflectors. Beams may be 
considered an oversimplification of the structural characteristics of segmented reflectors. 
However, the basic characteristics of sandwich construction are retained in the beam mod- 
els and beams could serve as a screening tools as effectively, and certainly as economically, 
as plate or shell-like models. Here efforts have been made to involve the important material 
properties explicitly so parametric studies can easily be made. Some approximations were 
necessary, but these have been justified and in no way do they compromise the results ob- 
tained. 

Several recommendations are in order. First, extending the analysis to include thermal ef- 
fects, such as would occur in the presence of a slight through-the-thickness temperature 
qradient, would be worthwhile. Such a gradient would cause unwanted curvature in a beam, 
and over time, the curvature may change. Second, it would be useful to extend the analysis 
to include the two-dimensional aspects of the reflector, namely its plate-like geometry. ” the 
change of focal length with time, for example, of a reflector is to be determined, such an 

analysis is necessary. 
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Fig. 1 - Geometry, nomenclature, loading, and assumed displacement field for the three- 

layer model. 
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Fig. 2 - Percent increase in tip deflection of three-layer beam due to time-dependent 
face-sheet extensional properties. 



Time (years) 

Fig. 3 - Percent increase in tip deflection of three-layer beam due to time-dependent core 

shear properties. 
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Fig. 4 - Geometry, nomenclature, loading, and assumed displacement field for the five 
layer model. 
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Fig. 5 - Percent increase in tip deflection of five-layer beam due to time-dependent ad 
hesive layer shear properties. 
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